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Abstract 

This paper is devoted to the statistics of the quantum eigenfunctions in an 
ensemble of finite disordered systems (metallic grains). We focus on mo- 
ments of inverse participation ratio. In the universal random matrix limit 
that corresponds to the infinite conductance of the grains, these moments are 
self- aver aging quantities. At large but finite conductance the moments do 
fluctuate due to the long range correlations in the eigenfunctions. We eval- 
uate the distributions of fluctuations at given conductance and geometry of 
the grains and express them through the spectrum of the diffusion operator 
in the grain. 
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Weakly disordered metallic grains make an excellent laboratory to study the phenomenon 
of Quantum Chaos (for a general discussion see e.g. Rcf. [|l|]). Provided electrons within a 
grain interact weakly, one can describe properties of this system through one particle quan- 
tum spectrum and corresponding eigenf unctions. The problem is reduced to a Shrodinger 
equation for a particle subject to a potential that consists of two components: a regular po- 
tential that confines electrons within the grain and some random potential due to disorder. 
Given the distribution of the random potential we get an ensemble of disordered grains and 
can consider various statistics of the spectra and eigenf unctions. 

Classical motion in a random potential is diffusive (provided the grain size R exceeds 
the mean free path) with a diffusion constant D. Ensembles of weakly disordered metallic 
grains can be characterized by the "dimensionless conductance" g determined as a ratio 
g = tn/tr of Heisenberg time = h/A and Thouless time of the diffusion through the 
grain tx = R 2 /D, where A is the mean energy level spacing. A grain would be called weakly 
disordered provided g 1. 

The spectral statistics for the ensembles of grains in the limit g —>■ oo are proven to 
coincide with those for the corresponding ensembles of random matrices [3J]. These statistics 
would be called universal. The field theoretical way of evaluating statistics in ensembles of 
disordered grains is based on supersymmetric cx-model 0. When g — > oo the O-dimensional 
ex-model can be used for straightforward evaluation of universal statistics of both spectra 
and eigenf unction ||. 

Finite g corrections to the universal properties of quantum systems recently attracted 
substantial interest. The smooth part of the spectral correlation function was evaluated in 
Ref. ||. The first order in 1/g correction to the spectral correlation function was evaluated 
by Kravtsov and Mirlin |7]]. Nonperturbative analysis of the problem || pointed out a 
qualitative change of the behavior at finite g - washing out the oscillations in the two-point 
correlation function, and, hence, smoothing of the singularities |9j in the quantum dynamics 
at times close to A remarkable feature of the spectral statistics at finite g is that 
both their smooth [f| and their oscillating parts f| can be presented through the spectral 
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determinant of the classical diffusion operator DV 2 . This paper is devoted to the connection 
of this spectral determinant to the long range correlations in the eigenfunctions. 



In this paper we concentrate on the moments [10 



I a (n) = V n - 1 J |V Q (r)| 2n rfr, (1) 

where ip a (r) is an eigenfunction of the system which corresponds to an eigenenergy e a , and 
V is the volume of the grain (for a d-dimensional cube V = L d ). The n = 2 moment, 
known as the inverse participation ratio, is related to the level- velocity distribution JO] or 
to Hubbard-like interaction of two particles on the same quantum sate. 

In the universal regime g — > oo each wave function is extended over the whole volume, 
however only very short range correlations persist: ip a {v\) and ip a (r 2 ) are correlated provided 
r = | i"i — r 2 | is not much bigger than the particle wavelength ||,|1^-[14|| . As a result, the 
integration in Eq.(|l]) provides self- averaging, and I a {n) do not fluctuate in the universal 
regime! They coincide with the moments b n of the Porter-Thomas distribution for the 
intensity fluctuations [^(r)! 2 - For unitary (it) and orthogonal (o) symmetries they are equal 
to I a {n) = b n , where 

&£•> = T(n + 1), 6f) = 2T(n + l/2)/r(l/2). (2) 

The most striking difference of finite g case from the universal situation is the existence of 
spatial correlations of wave function density even at r comparable with the size of the system. 
As a result I a (n) demonstrates finite fluctuations from state to state and from sample to 



sample [fLl|j . I a (n) are characterized by their distribution functions. Let us consider the 



distribution P n (u) of relative deviations of I a (n) from b n : 

u a (n) = I a (n)/b n - 1, P n (u) =<S(u-u a (n)) >, (3) 

where < ... > stands for ensemble averaging. 

As we show below the Laplace transforms P n (s) =< exp(—su a (n)) > of the distribution 
P n (u) for large but finite g can be written as 



-1/2 

pm = n 



1+ n 2 -n — — 



N «((»-»')^). W 

where /3 is different for T- invariant and not T- invariant systems: (3^ = 1, (3^ = 2 (we are 
not considering a simplectic case here), uj^ is the spectrum of the diffusion equation with 
the Neumann boundary conditions on the grain boundary B: 

DV 2 ^(r) = -u^r), V0| b = 0. (5) 

uj^ are not universal: this spectrum is determined by both g and the shape of the grain. 
(We take g = j3ui/(2A) as a definition of g, since this ratio is proportional to the dimen- 
sionless conductance for a rectangular grain). However, the ground state of the problem is 
spatially uniform, and corresponds to Uq = 0. All universal statistics become applicable to 
disordered grains when uo^ — > oo for all fi except fi = 0. In this limit, known as zero mode 
approximation, P n (u) = S(u) since Eq.(f|) gives P n (s) = const. Therefore I a (n) does not 
fluctuate in the universal regime. 

In Eq. (|4]) we introduced the function 

l/Z(z) = H[l-z/u> ll ] = z/Z(z), (6) 



where Z(z) is the dynamical Ruelle zeta function |T^], associated with the diffusion operator 
in Eq. As it is shown in Ref. |l6j the pair spectral correlation function for an ensemble 
of disordered grains also can be expressed through Z(z)-function. However the spectral 
statistics are determined by \Z(iz)\ 2 , while for the distributions P n (u) one has to evaluate 
Z(z) at real negative z or to determine both modulus and phase of Z(iz). An interesting 
feature of the distributions P n (u) is that, in contrast with other statistics of the quantum 
eigenf unctions, they are independent on the diffusion operator eigenfunctions M , and are 
determined solely by the spectrum of Eq. || 

Let us describe the main features of the distributions P n (u). It is clear that long scale 
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correlations increase the mean value < J(n) >= & n [l+ < u >] of J a (n) as compared with 
Eq.(g). According to Eq.(|) 



1 



(ri 2 — n) 

and for < « >= / uP n {u)du we obtain 
< u > 



P(w), w = u n /(n 2 - n), 



(7) 



< w > 



n A — n 



~~ 2 -j-Pn(s)\s=0 — -5— ^2 • 

n 2 -nds (Jn ±? lo^ 



In terms of the Zfzl-function < w > can be rewritten as 



s~> / \ /2A\ m tf m lnZ(z) 



(8) 



Behavior of P(u>) at small w can be evaluated by making a saddle-point approximation in 
the inverse Laplace transformation of Eq. 
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P( W ) = Z 1 / 2 (^)[vr^(^)]- 1 / 2 exp 



2A 



(9) 



provided z c (iu) determined by the equation Qi(z c ) = 2w is large, \z c \ ^> uj\. 

It follows from Eq.(f|) that the probability for w to be much bigger than < w > is 
exponentially small: 



P (w) =C Jg/(Aw) exp [-ngw], 



(10) 



^0,1 



-1/2 



I 2— >LJl • 



Consider a disordered two dimensional grain with a particle mean free path /. From 
Eq.(|) it follows 



0i(z) = -h 
5 



1 + 



tJi P 2 

0?i — -2 I 2 



< w >-- 



\n(R/l) 
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(11) 



Therefore in the weak localization regime [17] when all quantum states are extended < 
w ><^ 1. Eq. ( |TT|) for < w > is in agreement with the perturbation theory calculation [11 
For 1/jCii) <C< tu > Eq. (H) gives 



= I exp 



g(< w>-w)- |e 2s(<u,> - u,) 



where < w > is determined by Eg. (|TT|) . For gw < 1 



P(w) 



7 3/2 



exp 



g < w > 



71 



,2g<w> 



4gw 



According to Eq.(|10D InC = ng < w >, and at < w >< w <C 1 the distribution is 



P(w) = \Jg/(4w) exp [— Trg(w— < w >)]. 



(12) 



It should be mentioned that Eqs.([I|-|T2"D are valid only for u < 1. According to Eq.(^) 



this means that Eqs.(f|- |TJD describe the main body of the distribution P n {u) which never 
takes a gaussian form. 

When m>1, the distributions P n {u) are determined by both ui^ and 0^. This asymptotic 
of P n {u) can be found by the method of optimal fluctuation. For a spherical grain of arbitrary 
dimension d and radius R 



P n (u > 1) « exp[-a n gu 



l/(n-l)i 



(13) 



where # = ir/3D/(4Jl 2 A) and 



d 2 



-r 2 



1-2/ 



n 



n 



l/(n-l) 



7m(n-l) 2n-2 y y 2n-2 
This distribution is valid as long as u <C (a n R/l) n ~ l . At large w using the cx-model approach 
fails [0 . 

In order to derive Eq.(^) we consider mutual distribution of v\ = V^I^q, | 2 and 
v 2 = V\ip a (r2)\ 2 - densities of a quantum eigenfunction ip a {r). This distribution can be 
reconstructed through its moments M pq = V p+q < \ip a ( r i)\ 2p \i J a(' r 2)\ 2q >• The latter can 
be calculated for a disordered grain by analyzing moments of one-electron Green functions 
using the supersymmetric a-model technique in a way similar to how one-point moments 
M p o were calculated by Muzykantskii and Khmelnitskii 0] , and by Fal'ko and Efetov |f20 



The ratio of M pq and its universal (g — > oo) value M^ n ^ can be written [21j as the 



functional integral 



M 



i J V9(r)e 



-F+ P (e 1 -n)+g(9 2 -n) 



(14) 



where 9\ }2 = 9(ri t2 ) an d 5 = / T>9(j)e F . Note that while M^ 1 is determined by the zero 



mode of the diffusion operator, the ratio Eq. (H) does not include the integration over this 



mode. As a result the functional integral in Eq. (|bj|) is over all functions 9(r) that satisfy 
the condition / dr9 = 0, and 

^ = ^K «M = -(/4)- < 15 > 

Eq.(^4|) enables us to 

express the two-points mutual distribution function P(V|-?/> Q (ri)| 2 , V\ip a (r 2 )\ 2 ) at finite g 
through the universal one P un { v ii v 2) 

P( Vl ,v 2 ) = 1 / V9(v)e- F+2n - e ^ 

xC^e^',^). (16) 

The universal two-point distribution functions P un {vi,v 2 ) for unitary and orthogonal 
symmetries were determined earlier |j ,(Li|]: 

V1+V2 



, \ ex p 



2(1-/ 2 ) 



/n 1 1 _ f2 

f^f>HV2 



t>l+t>2 

cosh 



2^^/(1 - p) Vl v 2 V 1 -/ 2 /' 

where / is the Friedel function of the distance r = |ri — r 2 \ 

f(r) = V{d/2){2/kr) d ' 2 - 1 J d/2 -x(kr)e' r ^ l \ 

Above J p (x) and Iq(x) are Bessel and modified Bessel functions, respectively, k is the wave 
number. Note that f(r) — > 0, when kr 3> 1, and /(0) = 1. 

A usual way to calculate the functional integral like Eq. ( |16"D is to present 9{r) as a sum 
S/i^^( r ) (there is no contribution with ji = since J dr0 = 0) over the eigenfunctions of 
the problem Eq.(^) thus reducing the calculation of the functional integral to a sequence of 
definite integrals over 9^. 

Expansion of the exponent in Eq. (|T6| ) up to the second order in 9^ leads to a gaussian 
integral that can be evaluated explicitly 



P(vi,v 2 ) 



OO /"OO 



dsids 2 Pun(v 1 s 1 ,v 2 s 2 ) 



o Jo 2ti 



n n n 22 - n 2 



.12 



x exp 



E 



(vTTn~lns2 + «Vn 22 lnsi^ 



K=±1 niin 22 + Kni 2V /nnii 22 

where Iljj is the Green function of Eq.(|5|) 

For ri — > r 2 Eq. (|17D reproduces the result for one-point fluctuations pO 



(17) 



(18) 



From Eg . (|16D we can guess the probability density of a given realization ip(r) of an 
eigenfunction in the form of a functional integral. Instead of the gaussian distribution that 



is valid in the universal limit 13,14 we obtain 



P[$\Thl) = Vi\>- J V0(r)e~ F (/T 1 det K) P 



x exp 



2 



(19) 



The matrix elements of the operator K in the coordinate representation fT 12 =< r 1 |i^|r 2 > 
equal to 

0i + 2 



K12 = /(|ri - r 2 |) exp 



(20) 



The distribution determined by Eqs.(19,20) enables straightforward calculation of cor- 
relation functions for wave function intensity at many points separated by distance larger 
than the mean free path: 

< n^(r,)| 2 "< >= 4/ D0(T)e- F Hb ni K%- (21) 

i ' ' i 

For two points correlation (i = 1, 2) Eq.(21) directly follows from Eq. ([14]) . The many points 



correlation functions also can be shown to coincide with Eq. (21) 
Eq.(|19|) enables us to present P n (s) in the form 



Puis) 



V0{v) ^ F 
— ^ — e ex P 



l)dr 



(22) 



Using in Eq. ( p2j) the harmonic approximation that is valid for small fluctuations u <C 1, 
we obtain Eq. 

At m > 1 we can apply the method of optimal fluctuation to Eq. ( fffi) . For a spherical 
dot of the radius the saddle-point equations can be written as 

t d_1 dtj, (23) 

where t = r/R, a dot stands fot t-derivatives, and 9 c (t) obeys the equation 

e c +^e c = -~ $(0)=p(e" e -n e e ), (24) 

supplied by the condition 9 C (0) = 9 c (l) = 0. p in Eq. (|24D should be determined from the 
self-consistency condition 

1 — n V— sj Jo 

Let us note that the potential has a minimum at 6* = 0. At u ^> 1 or \s\ <C g 
the essential contribution to Eq. (^) is given by the trajectories of 9 c (t), that starts with 
9 C {0) ~ — ^2 l n l^/ 5 ! ~~ > — 00 an d ends somewhere near the minimum of potential They 
correspond to the optimal wave functions that in the center of dot is small like \s/ g\ 2 ^ n ~ 2 ^ 
and increases when approaching the boundary. Such a trajectory leads to the asymptotic 
Eq.(|i~3"D for the distribution function P n (u). 

Fal'ko and Efetov |22[] found that there are long range spatial correlations of eigenfunc- 
tions even at g — oo, provided in a crossover from orthogonal to unitary symmetries, rather 
than any of the pure symmetry cases. Here we have not discussed the crossover regime. 
Calculations similar to ours for the distributions of moments Eq.(|14"D and the two point 
distribution function Eq.flTTD in the crossover regime can be carried out rather straightfor- 
wardly. Details of this calculation will be published elsewhere |[21|| . The statistics for the 
crossover regime are determined by so-called "Cooperon gap" u c ~ ($/$o) 2 ^i ; where $ is 
the magnetic flux through the system, and $ is the flux quantum. If u>x — > oo and u c 



exp 



gd 



91 + 



n — 1 
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is finite, these results are in agreement with Ref. |22j. In the crossover regime the diffu- 
son eigenf unctions are not orthogonal to cooperon ones. Therefore the distributions of the 
moments Eq.(|14]) acquire dependence on the eigenfunctions. 

To summarize, we have calculated the distributions for generalized moments of inverse 
participation ratio Eq.(l) in an ensemble of disordered metallic grains with a given dimen- 
sionless conductance g. In the universal limit g — > oo these moments do not fluctuate due to 
self- averaging and thus have definite values Eq.(2). The fluctuations appear only at finite g 
together with long-range correlations in densities of wave functions. Indeed, in the universal 
limit, when only short-range correlations persist, the fluctuations vanish when the grain 
volume (the rank of random matrices) tends to infinity. Contrarily, at finite g they appear 
small (as 1/g) but long range correlations (see Eq. fll8|) ) controlled entirely by the eigen- 
functions of the diffusion operator Eq.(5). These correlations give rise to the fluctuations of 
the moments. 

It is amazing however that at g ^> 1 and for n <C ^fg the main body of the distribution 
functions P n (u) (Eqs. (|T|,|3|, |l0| - |T2|) ) can be expressed through the spectrum of the diffusion 
operator and do not depend on the eigenfunctions This fact suggests generalization 
of these statistics of eigenfunctions from disordered to generic chaotic systems, by making 
use of the Ruelle (dynamical) zeta-function (Eqs. f|||), in a way similar to what was done 
in Ref. |L6f for spectral statistics. Substitution of the spectral determinant of the diffusion 
operator by the Ruelle zeta-function for spectral statistics in generic case was supported by 



the calculations in the framework of a nonlinear cr-model |23]. The analysis of Bogomolny 
and Keating |24| based on the periodic orbit theory led to similar but different results. We 
hope that further analytical development of results on the eigenfunction statistics together 
with numerical evaluations will clarify the relation between quantum statistics and classical 
dynamics including the phenomenon of scarring |[25|| . 

The authors are grateful to O. Agam, I. Aleiner, A. Andreev, and K. Efetov for useful 
discussions and thank P. Fulde for his interest to the work and support. 
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